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PREFACE 


The tumbling and purface roughness effects on the trajectory 
of entry tektite are studied in both free molecular and continuum 
flows. It was concluded thatj while surface roughness has neg- 
ligible effect on trajectory^ the tumbling may play an important 
role in tektite trajectory and the consequent ablation, provided 
the body shape is different from a sphere . A shape factor B was 
proposed to measure the shape irregularity and was found to be 
a good parameter for correlations between body shape and tumbling 
effects . 
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Nomenclature 


area 

deceleration 

thermal accommodation coefficient 

reference cross-sectional area 

shape factor 

drag coefficient 

skin friction coefficient 

diameter of equivalent sphere 

cylinder diameter 

axial drag force 

total drag force 

transverse drag force 

coordinate vector in Pig. 1 

coordinate vector in Pig. 1 

coordinate vector in Pig. 1 
roughness height 
Knuds en number 

unit vector along cylinder axis 

clinder length 

length scale 

Mach number 

pressure 

dynamic pressure 
Reynolds number 
nose radius 


Nomenclature (Cont'd) 


total surface area 
temperature or total time 
time 

velocity 
total volume 
weight 

angle of attack 
ballistic coefficient 
flight path angle 
shock stand-off distance 
boundary layer thickness 
density ratio across the shock 
speed ratio 

angle between V and W (Pig. 1) 

azimuth angle of cylinder axisj or angle between 
surface normal and velocity 

mean free path 

viscosity 

density ■ 

angular velocity 
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Subscripts 
A axial 

C continuum 

e entry condition^ or edge of boundary layer 

ed end-on position 

PM free molecular flow 

/I roughness 

s stagnation, or side wall 

sd side -on position 

sp equivalent sphere 

sm smooth 

T transverse 

w wall 

z conditions behind shock 

00 free stream 
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INTRODUCTION 


Most works of tektite ablation studies^ ^ ^ were under the 

assumption that the entry tektite is of spherical shape and smooth 
surface. The numerical results of those calculations indicated 
that ablation takes place mainly by melt flow rather than by vap- 
orization. {1“2) This is disturbing in comparison with the observa- 
tions on the recovered tektite, since melt flow has not been definit- 
ely found on tektite except for autstrllte and a few Javanites . 

The effects of surface roughness on boundary layer transition, 
heating and ablation have been studied in Reference 2 >, which con- 
cluded that the surface roughness on most recovered tektite was 
formed before earth entry. Consequently, the roughness may trip 
the boundary layers to turbulence and enhance vaporization. Num- 
erical results in Reference 3 indicated that the effect of surface 
roughness on ablation appeared to be an important factor to ex- 
plain the discrepancy of vaporization vs melt flow. The effect 
of surface roughness on drag has not been discussed in Reference 3. 

Most of the entry tektites are expected to be nonspherical in 
shape and to carry a not negligible rotational momentum. In 
addition, irregular or asymmetrical bodies, in which the resultant 
of aerodynamic pressure doesn't pass through the center of mass, 
will experience the action of a variable couple influencing rota- 
tion. Thus, one way or the other, the entry tektite is expected 
to be tumbling in some fashion. Since the drag coefficient is 
highly dependent on the angle of attack, the effective ballistic 
coefficient (or the trajectory) and the total ablation may be sub- 
stantially different from the value of corres poinding spherical 
bodies . In this paper, analysis were conducted under the assump- 
tion that the frequency of tumble is large compared to the change 
in altitude during any cycle; thus, the average drag coefficient may 
be associated with the fixed environment of a single altitude. 

The correction of drag coefficient will change the trajectory of 
entry tektite, affecting the corresponding albatlon behavior. 

During earth entry, tektite experienced various regions of 
atmosphere, from free molecular flow, transition flow to continuum 
flow, depending on the corresponding density level. The division 
of gas dynamics into various regions are usually conducted based 
on characteristic reanges of values of an appropriate Knudson 
number, defined by a dimensionless ratio X /L, where x denotes 
molecular mean free path and L is a length scale. There are 
various ways in defining X and L; e.g. the mean free path 
may be measured at freestream or behind the shock, while the length 
scale L may be referred to as the nose radius boundary layer 

thickness b , or the shock stand off distance A, depending on 

the parameters and the environments. 
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Plow in the various transition regimes between free -molecular 
and continuum is extremely complicated^ and no satisfactory 
theoretical solution has yet been offered. No detailed develop- 
ment of drag coefficients or skin friction will be presented in 
this report for transitional flows. Discussions will be concentra- 
ted on the limits of free-molecular and continuum flows. Engineer- 
ing methods will be presented for extrapolating drag coefficients 
and skin frictions between these two well defined regimes . 


In free molecular region, the mean free path is everywhere 
very much greater than a characteristic body dimension, and the 
distortion of the free stream velocity distribution due to the 
presence of the body is negligible. However, the actual transfer 
processes are tied up with the fact that the molecules are remit- 
ted or reflected in some manner j thus, a thermal accommodation 
coefficient is introduced to measure the extent to which the mean 
energy of the remitted molecules is accommodated toward an energy 
corresponding to the temperature of the wall. The thermal accom- 
modation coefficient is defined by the equation 


E; - Ea 
^ e EJL - E b- 


( 1 ) 


where EX is the incident energy per unit surface area per second, 

Ea is the reflected or re-emitted energy carried away by the 
molecules as they leave the body, and Ey is the energy that 
the re -emitted or reflected stream would have if all the incident 
molecules were re -emitted with a Maxwellian velocity distributions 
corresponding, to the surface temperature T , . Tests of Wledmann 
and Trumplerl^) indicate that under statlc’^condltlons, the thermal 
accommodation coefficient for air on various typical engineering 
surfaces lies between 0.87 and O.97. Thus, it is reasonable to use 
unity accommodation coefficient for most practical problems . 

The other important parameter in rarefied gases is the so 
called speed ratio, f , which is defined as the ratio of the mass 
velocity to the most probable random speed in the freestream. 

It was estimated in Reference (5) that the average of air molecules 
at high altitudes to be about 4000-5000 ft ./sec.; an order of magni- 
tude less than the typical speed of tektlte entry. Under the assump- 
tion of unity accommodation coefficient that the air molecules are 
reflected at roughly the same temperature as the tektite, the re- 
flected particle should have velocities at the neighborhood of 
3000 to 4000 ft ./sec. Since the molecular velocity before and 
after the collision with tektlte is much less than tektite velocity, 
the molecular velocity incident to tektite can be approximated 
by the tektite velocity itself and the rebounclng velocity of air 
molecules can be approximated by zero. Thus, the approximation 
of infinite speed ratio, will be used throughout the analysis 

of free molecular flow . 
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The obher limit of the flow regimea la the well-studied con- 
tinuum flows at sufficiently high Reynolds numbers . The viscous 
effects in this region may be taken into account using classical 
boundary layer theory, with suitable modification for pressure 
interaction if appropriate. Pressure drags on bodies traveling 
with hypersonic speed can be estimated with Newtonian approximation, 
while surface roughness effects can be calculated from classic 
boundary layer theory. With the free molecular and continuum flow 
results calculated, the quantities in transitional flows can be 
interpolated between these two well defined regions . 

Tektite analysis contain a source of error depending on the 
various shape of the tektite, except when the latter is nearly 
spherical such as australlte. However, australlte composed only 
a small portion of the recovered tektite, and indeed most of the 
tektite are irregular in shape and do not lend themselves to purely 
analytical analysis. Certain plausible assumptions can be made in 
this respect for order of magnitude calculations. To measure the 
geometry of an entering tektite, a shape factor b(o) is defined 
as the actual surface area divided by the surface area of a sphere 
of equivalent volume, 

B = S/r (2) 

where S is the tektite surface area and D is the diameter of the 
equivalent sphere, 

D = (6 (3) 

with Vp denoting the volume of tektite . 

Parametric studies for tektite shape effects on trajectory 
and ablation will be carried out in terms of shape factor B with 
circular cylinder used as basic shape for numerical results . It 
is hoped that simple analysis can bring cut the order of magnitude 
effect of body shape on tektite entry behavior. 

The surface area and volume of a circular cylinder of diameter 
d and length i are respectively as follows: 

S = TT/2 d^ + 7T dfi 
Vp= *^^4 di 

Substitution of Eq, (4) into (2) and (3) yields the diameter of 
equivalent sphere D and the shape factor B, 

D ~ 1.1447 i 

B = 0,3816 + 0.7631 (i/d)^/^ 
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Dependence of B on d/j2 is plotted in Figure 1. It is indica- 
ted that B tends to increase in both long cylinder (£/d > l) and 
short cylinder (£/d < 1) with a minimum at j? « d. Both long cylinder 
and short cylinder results will be presented in this report, and 
the Inference of these results to tektlte study will depend on the 
actual shape of tektlte . 
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TUMBLING AT FREE MOLECULAR FLOWS 


The trajectory of a body entering the atmosphere from high 
altitude depends greatly on the entry velocity Ye, the entry flight 
path angle and the ballistic coefficient ^ which is defined as 

w/CjjA (6) 

where W is the weight of the body^ and and A are the drag coef- 
ficient and the corresponding reference area respectively. While 
Ve and specify the initial conditions of earth entj'y^ the bal- 

llstic coefficient is a function of Mach number, Reynolds number, 
and geometry of the entry object . Thus, the key factor for trajectory 
calculation is the determination of 

Under the approximations of unity accommodation coefficient A„ 

and infinite speed ratio , the drag on tektite at free molecular 
flow is actually equal to the change in momentum of all the air 
molecules striking the surface per unit time j i.e. 

= Po. ^ ( 7 ) 

where and A denote the corresponding atmosphere density, 

tektite velocity and projected front area respectively. Thus, Cj^ 
for free molecular flow can be well approximated by 2, based on 
the projected front area. This leads to the following expressions 
of ballistic coefficient 

= W/(2 d) (8) 

for cylinder at side-on positltion {cylinder axis_L flow) 

^ed SW/Ttd^ (9) 

for cylinder at end-on position (cylinder axis II flow) and 

/3gp = W/(2.0583 d^/3 £^/3) (10) 

for equivalent sphere, where W denotes weight, and I and d are 
the length and the diameter of the cylinder respectively. 

The same argument leads to the corresponding ballistic coef- 
ficient 

0= W/(2d£ sin a + Jj/2 d^ cos a ) (11) 

For a tumbling body, the average ballistic coefficient is defined 
as being w/C^A of an object which has the same deceleration 


_ 5 _ 


i 


as the average decelerations: 


a = 1/T J a ( t ) dt 


( 12 ) 


where a and t denote deceleration and time respectively, and 

is the dynamic pressure which is assumed not vary greatly during T. 
By definition of average W/Cj^A 


® = 9<x./ ^ 

and, therefore. 


_ T Jo 


(13) 


( 14 ) 


If the motion of^ a cylinder is assumed tumbling in the trajectory 
plane (i.e. w V = 0 / and the variation of angular velocity ui 
is small during T, thus. 


= wt 


(15) 


and the average ballistic coefficient in terms of is given by. 




substitution of Eq. (ll) into eq. (l6) yields 
W/(VlTid + d^) 

However, since tektites are probably put out at random spin 
orientation, it will in general be true that 

(j • 


(16) 


(17) 


(18) 


so that the actual will be different from that given by eq. (17). 

Figure 2 ^hows a {JL coordinate system with velocity 

vector in the k direction, _yie ^vector at an angle 9 with the 
velocity vector and in the k - J plane. The angle_oc is the angle 
between the cylinder axis and the velocity vector V, an^the 
angle 0 = iot is the azimuth angle of the cylinder axis,X . Prom 
Figure 2 it is obvious that 


= k cos e + 
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( 20 ) 
( 21 ) 

Substituting (21) inco (11), and (l6) we obtain the approximate 
solution 


i ^ x-g- ) cos + Zsln j?f 
coc oC = k: • £:= ain 0 sin 0 


' ( 


^ (0) = W/ (2 - -I sin^0 - 3/32 sin^0) + d^ sin0 


( 22 ) 


where a slight adjustment has been made to some of the sine terms 
to allow ^(77/2) matching eq. (17). Thus, ar each value of 0 = 0 * 
(0*^ ©■‘^ %'^ 2 ) an offect ballistic coefficient >B(0 q) can be ° 

determined which represents the average taken over a cycle of 
tumbling in the plane defined by © . As indicated earlier, no 
attempt is made to predict specific orientation of the plane of 
tumble in the present paper. Instead, an average value of the 
effective ballistic coefficient is defined for a random orienta- 
tion of the Plane of tumbling by the following equation. 

|=(E = 1 ! CS) '=1 S (23) 

Where^ is the probability density. Under the assumption that a> 
is ran(5omly distributed in space, it can be easily shown that the 
distribution function for 0 is given by 


P (©<0 q) = 1 - cos 0 q (24) 

and the probability density is given by, 

-^p 

Substitution of Eq. (25) and (22) into (23) yields, 

g =W/(l. 6l66 jgd + 0.7854 d^) _ (26) 

• * 

Equations (8), (9), (17) and (26) represent the ballistic 

coefficients of a cylinder shape at various position and tumbling 
patterns. Their relative values with respect to that of an 
equivalent spheh* are summarized in Table 1. Combination of 
Table 1 and eq. (5) yields the numerical results plotted in 
Figure 3, which shows the ratio of ballistic coefficient to that 
of an equivalent sphere as a function of the shape factor B. 


Since the shape of long cylinder differs very much from that 
of the corresponding short cylinder of the same shape factor B, 
the relative ballistic coefficient of an irregularly shaped tektite 
is expected to be somewhere between or close to those of the long and 
short cylinders. Both results of tumbling in pitch plane and of 
random tumbling are included in Figure 3 to cover the range of 
various styles of tumbling. Thus, once the volume and the surface 
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area of a tektlte shape are measured, the range of ballistic 
coefficient can be estimated from the cross-hatched area in 
Figure 3 . Since the shape of tektite can have a shape factor B 
of 2~3j the drag acting on it can easily be 2—3 times larger 
than that of the equivalent sphere. 
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TUMBLING IN CONTINUUM AND TRANSITIONAL PLOWS 


As the tektite approaching continuum flow region, the mean 
free path becomes negligibly small compared to the dimension of 
tektite and the viscous effects are limited to a thin layer over 
the body surface (boundary layer). Because of the hypersonic 
speed, a modified Newtonian pressure distribution appears to be 
appropriate for mose practical problems. The local pressure is 
given by 

f = ^00 + - J'oo) 

where P is pressure, ^ is the angle between the surface normal and 
the velocity and subscripts oo and s denobe the conditions at free- 
stream and at stagnation point respectively. The continuity and 
momentum equations across a normal shock are 

P 00 ^00 "" ^2 ^2. 


^00 + 


? V 

GD oo 




2 

2 ''2. 


V, 


Where pand V are density and velocity respectively, and subscript 2 
denotes condition immediately after the normal shock. Since 
incompressible flow relations closely approximate actual conditions 
in the stagnation region behind the shock wave, the stagnation 
pressure can be given by 







Thus, eq. (27) becomes 


P + =100 (2 - Pj p.^) 


For Mach number greater than 10, P_, q^ and P„ ^ — .q_^ (2 - 

00 00 S — 00 ' oo' d. 

and the following simplified relations have been frequently used 
for practical purpose: 


P ^q^ (2 - €) cos^ ^ 


( 28 ) 


^ P_ cos^ ^ 

where £ denotes 9^/ p,^. 

The drag coefficient is given by 

Cp = y? ^ ^^/^oo 


(29) 


which is integrated over the surface of the body exposed to the 
freestream. For the case of a cylinder in side-on position, 

Eq. (29) yields 

Od = 2/3 (2 (30) 

which lead the side-on ballistic coefficient 

id (a -g )| (31) 

In addition, Eq, (29) leads to the drag coefficient and ballistic 
coefficient of the equivalent sphere as follows: 


Op =• 1/a (a -£ ) 


,V3 


^2/3j 


(32) 


^’sp “ fo-sine (a - e) d‘ 

As to the end-on position, there would be a uniform pressure P 
over the face of the cylinder according to Newtonian theory. ® 
However, Stoney and Swansom?) found that the average pressure on 
the face of a cylinder O.909 times . stagnation pressure. Using 
this value in Eq. (29) and integrating, we have, for a cylinder 
in end-on position. 


Op = 0.909 (a -e) 


(a - e ) a®] 


(33) 


For a cylinder at angle of attack, the total drag force can be 
approximated by the vector sum of the axial and transverse com- 
ponents 

Pp ■■= P^ sin Fj^ cos d 


or 


but 


and 


^Doc ^00 "A 




sinc»+ V « 


^ ^ foo ^00 


" H cos^ 

^A J CO A 00 


(3^n 

(35) 

(36) 


Combination of eqs . (3^)^ (35)3 (36), and the definition of bal- 
listic coefficient yields. 


^ (cL) . W/ 


sinocid + 0.71^ cos cx. d) (2 - £) 
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Vi ‘i* '• "nil" 


The average ballistic coefficient for a cylinder tumbling in the 
trajectory plane is given by substituting Eq, (38) into ^16), 

>0 = W/ jCo. 2829 id + 0.3030 d®) (2 - (38) 

For cylinder tumbling in plane at an angle ( *^/2 - 0) with the 
trajectory plane, the average ballistic coefficient can be 
calculated by the subst'’ tution of Eqs. (21) and (38) into (16). 

An approximate solution was obtained, 

(0) = w/ [^(0.6667-0.48 sin2 0 + 0.1 sin'^ 0 + 0.303 j sin^ G) (2-£) 

Where a slight adjustment has been made to some of the sine terms 
to allow^( “^/ 2 ) matching Eq. (38). The average balllstl * coef- 
ficient of a cylinder tumbling in random is given by subsDitutlon 
of Eq. (3? ) into (23), 

^ = W/ [(0.4 id + 0.1785 d^) (2 -£)] ( 40 ) 

Table 1 summarizes the relative ballistic coefficients of a 
cylinder traveling in continuum flow at various positions and 
tumbling patterns. Combination of Table 1 and Eq. (5) yields the 
variation of those relative ballistic coefficients against shape 
factor B as plotted in Figure 4. 

As for drag coefficient in transition regions ^.between free 
molecular and continuum flows, Bloxson and Rhodes(°' conducted 
a series of drag coefficient experiments in a hypersonic wind 
tunnel on bodies of various shapes with Knudsen number range 
0.0001 - 0.34. The Knudsen number there was defined as 1/3 mean 
free path behind shock divided by the shock stand-off distance, i.e. 

K„=A/ 3 A ( 41 ) 

The drag coefficients so obtained on spheres are plotted 
in Figure 5 as a function of Knudsen number, which shows that drag 
coefficients can be presented as an exponential function of Knudsen 
number in transition flow. 


The most remarkable point about Figure 5 is that all of the 
shapes considered have the same variation with K^ in the transi- 
tion region and the same drag coefficient in the free molecular 
region. The variation seems to hold for any shape and gas, having 
only the requirement that the flow is over Mach number 4 . An 
empirical formula was derived based on the data of spherical 
bodies in Reference 8; i.e. 



^D,C ^D, PM . 

73 f- \ 



1 52 

K^ * - 0.0504 

j^l.52 ^ 0.0504 


( 42 ) 





I 

where Gj^, ^ and denote the drag coefficienta in free- * J 

molecular and continuum regions respectively. 'I 



ROUGHNESS EFFECT 


Surface roughness is known to promote laminar- turbulent transi- 
tion and to increase surface drag and heating by a significant 
amount . The surface drag on rough surfaces depends on the number 
of roughness elements per unit area, on their shapes and heights, 
and on the way in v;hich they are distributed over the surface; 
since the main contribution of surface drag is from the pressure 
force acting on the roughness elements. The number of parameters 
required for roughness description is extraordinarily large and no 
analysis have been developed to include all the parameters required. 
Most authors measured the roughness effect on surface drag by a 
single parameter of effective roughness defined as the rough- 

ness size of closely packed sand roughness which gives the same 
value of skin friction coefficient as the actual roughness. 

In this way, the difficulty has been shifted to the deter- 
mination of effective roughness from the measurement of actual 
roughness . However, for order of magnitude analysis such as 
tektlte studies, the average height of the reasonably spaced rough- 
ness elements can be used as effective roughness for practical 
calculations. As discussed in Reference 3j an effective roughness, 
Ke£.f j equal to 40 mils is assumed for the tektlte studies in this 

report . 

The exact way to assess the surface roughness effect on 
trajectory is to calculate the drag coefficient of a roughened 
tektlte which is tumbling during entry. However, the coupling 
phenomena of roughness and tumbling is a very complicated problem 
because of the cross flows Involved in an Inclined cylinder, which 
is beyond the scope of this report. The approach here is to 
estimate the roughness effect according to the most noticeable 
case, l.e., the case that roughness has maximum effect. For a 
long cylinder (d/i<l), the roughnesp effect is maximum as the 
cylinder is at end-on position, and the numerical results of 
roughness effect will be carried out In this position. 

Roughness effects on drag are closely related to the behavior 
of skin friction drag coefficient through various regions of 
earth entry. Reference 9 presented an empirical relation of the skin 
friction coefficients in the form of (C™ M) as a function of the 
interaction parameter m/ VHe", which is ^defined for flows with 
significant viscous effect. In such a case, the significant 
characteristic dimension of the flow field is the boundary layer 
thickness 5 rather than a dimension L typical of the body itself, 
and the corresponding Knuds on number based on the length scale 5 
is actually proportional to the interaction parameter fi/ . 
Hoerner(9j suggested that for skin friction coefficient simple 
gas dynamics for continuum flows are applicable to m/V^< 10' j 
and a transitional phase between 10"2 and lCH-1, above which free 
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molecular flow 1b finally established. Figure 7 presented both 
laminar and turbulent skin friction coefficient in terms of 
C^M vs m/ %^Re. Very good correlation is found between free molecular 
and continuum results . Naturally, turbulent friction should not 
be expected to reduce below the level as found for laminar bound- 
ary layer flows. The experimental results as plotted in Figure 6 
suggest that the laminar function represents the limiting condition 
for the turbulent function. It appears that a low Reynolds 
number the distinction between laminar and turbulent boundary 
layer flow becomes lost. It is suggested\°/ that the vortlcity 
originating across a boundary layer induced shock wave behaves like 
turbulent, thus possibly rendering a "laminar" layer effectively 
turbulent . For surfaces containing roughness of 40 mils such as 
observed on many recovered tektite, the bou ndar y layers will be 
tripped to turbulence somewhere between m/ =slO"l and 1 
according to the numerical results in Reference 3; and turbulence 
exists in boundary layers in most of the transitional region. 

Thus, it is plausible to use the conventional analysis of rough- 
ness effect in turbulent flow through the whole trajectory of 
tektite entry. 


The effect of roughness on turbulent skin friction has re- 
ceived considerable investigation in low speed flow as summarized 
by Clauser.^^^^ The roughness skin friction in high speed flow was 
studied by Goddard, v/ho shov;ed that the effect of compressiblility 
can count as a reduction wall density as Mach number increased. 

In addition, t,ie skin friction increase due to roughness correlated 
with the roughness Reynolds number 

Independent of Mach number. The data of Reference 11, which were 
limited to adiabatic flow conditions, yield a close approximation 
for roughness effect on skin friction: 




(43) 


.( 12 ) 


More recently Young' ' found that as the wall temperature was 
reduced, the skin friction dropped below the Goddard value for 
adiabatic walls . Reasoning that the effect of heat transfer was 
to produce a density gradient near the wall such that the effective 
density at the roughness surface was less than the wall value. 

Young developed a reference temperature method that correlated his 
results. 


(Te/Tr) + 0.635 (Te/Tw) (44) 

where C^. is the incompressible value of C- including roughness 
effect. ^ 


The relation between roughness skin friction coefficient and 
its adiabatic wall value can be derived directly from Eq. (44) 

i ,e . 



V 
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= 0.365 (Tw/Pr) + 0.635 


(4-5) 


V°f,av 

Combination of Eqs. (45") and. (43) yields the relation be- 
tween rough wall skin friction ihcluding heat transfer to that of 
stnmooth wall under adiabatic conditionsi^S / 

C^/Cf = [ 0.365 (Tw/T./l) + 0 . 635 ] log 10 R6 j^* (46) 


Thus, once the reynolds number around the roughness elements, 
the roughness size K and the temperature ratio Tw/Tr are given, 
the affect of surface roughness on skin friction can be calculated. 
The skin friction coefficient of smooth wall can be estimated 

from the correlations in Figure 6. ® 

The total drag for a long cylinder traveling in an end-on 
position can be expressed as follows: 


Ft 


<^00 ^e 


(C 


f q 


00 




m 


where C^ term denotes the friction drag acting on the sidewall of 

the cylinder as indicated in Figure 6 and 0^ term indicates the 

pressure drag acting on the end face of the cylinder. Usually the 
friction drag for a smooth cylinder is at least three orders of 
magnitude less than the pressure drag, providing the area rateo 

is of order one. This explains why only pressure drags are 

considered in tumbling studies. For rough surfaces, the friction 
drag can be larger than that of smooth wall, and its effect can 
be expressed by the relative drag of rough to smooth long cylinders 
as follows; 



r-sm 




P. 


D, sm 




r 


s 


+ 



(48) 


where can be estimated from Eqs. (9)j (33)# and (42) for free 
molecular, continuum, and transitional flows respectively, 0^ 

from Figure 6, from Eq. (46)# %/%q hypersonic results 

such as in Reference 3, and A^/a„ from geometry considerations. 
The area ratio A_/A_ can be related to the shape factor B through 

S' © 
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Eq. (5). 

Figure 7 shows the relative drag of rough long cylinder to 
smooth one at end-on position. In which drag Is presented as a 
function of shape factor B. The Mach number here Is a constant 
value of 30 > which Is the average value In the actively ablating 
region of tektlte entry. (3) 

The Reynolds number here, based on the equivalent diameter 
of one inch, is lo5, which Is the upper limit that an entry 
tektlte may experience according to the numerical results In 
Reference 3. The ratio of roughness siae to the equivalent diameter, 
K/D, Is selected to be 0.04; since the average roughness on tektlte 
is about 40 mils and the average diameter of tektlte is probably 
at the order of one inch. Figure 9 indicates that the drag In- 
crease due to roughness at Rep - IQP is less than 10^ for usual 

tektlte geometry (say KB<37* Since roughness effect on drag 
increases with Reynolds number, the drag increase Indicated in 
Figure 7 gives the upper limit of roughness effect on drag for 
tektlte entry. The Reynolds number of an entry tektlte usually is 
much less than 105; thus, it is plausible to assume that the 
surface roughness effect on drag can be neglected for tektlte 
entry studies , 
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TRAJECTORY AND ABLATION 


The trajectory of a vehicle entering the atmosphere from great 
altitude depends greatly on the entry velocity the flight path 

angle Y' , and the ballistic coefficient's . In general, both the 
flight path angle and the ballistic coefficient vary due to the 
presence of gravity force and to the change of drag coefficient 
during flight. Hence, it is clear that exact solution of the traject- 
ory is formidable. However, in the trajectory region of active 
ablation, the velocity of the entry tektite is so high that its 
total drag is much larger than the gravity force and that the drag 
coefficient is nearly constant because of the corresponding high 
Mach number environment Thus, for a simplified trajectory 
analysis like this, it is appropriate to assume that both S' and;9 
are constant. In addition, the atmosphere density^ , which 
is a function of altitude only, is given by the rel^ion, 

/oo^o = ) (49) 

where P and X are constant equal to O.O9I lbs. /ft. ^ and 23120 ft. 
respectively according to the curve fit for high altitude atmosphere 
in Figure 8. Based on these assumptions the relations, Allen and 
Eggersl^^^ derived the following relation 


V = Vq exp 


-fco ''■/ (2 ) 


(SO) 


which can be used to estimate the vehicle velocity as a function 
of altitude once the parameters of the entry conditions, such as 
B given. 


Numerical results in Reference 3 indicate that ablation rate 
was practically zero if velocity reduced to or below 10,000 ft/sec . 
Thus, one can define a consolidating altitude at which the velocity 
Just reaches 10,000 ft/sec and no significant ablation there- 
after. Combination of Eqs . (45), (^0), and (5) yields the 
numerical results plotted in Figure 9 which show the consolidating 
altitude as a function of shape factor B. The ballistic coefficient 
foimiula used in Figure 9 is based on that of random tumbling in 
free molecular flow with the equivalent diameter of one inch. 

It is clearly indicated that the shape factor B has significant 
effect of tumbling on trajectory. In addition, the geometry style 
(long or short cylinder) has only slight effect on the consolidating 
altitude, indicating that shape factor B is a fairly good para- 
meter to describe Iri’egular shapes. Figure 10 shows the similar 
conclusion for the trajectory calculations based on the randomly 
tumbling ballistic coefficient in continuum flows. 
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‘ r an order of magnitude estimation, the total recession of 
an ei y body can be measured by the total heat transfer to the 
body, 'fhe stagnation point heat transfer rate to a sphere of 
radius R 
equation 

* o 

where the unit of heat transfer rate q^ is in BTU/ft -sec, radius 

Rj^ in feet, and velocity V in ft/sec , Since the most obvious 

effect of tumbling on tektite behavior is velocity reduction due 
to higher drag, it is plausible to measure the ablation rate by 
the substitution of reduced velocity into Eq. (5/). The relative 
total ablation of velocity reduced tektite with respect to that 
of equivalent sphere can be measured by the total heat transfer 
to the stagnation point, i.e., 

A^tumb ^^tumb^ 

whereA S is total recession, subscript tumb and o denote tumbling 
correction and no tumbling respectively. Combination of Eqs . (51), 
{Si), and (5) yields the numerical results plotted in Figure 10, 
where the relative total recession is presented as a function of 
shape factor B. The ballistic coefficient formula used in Figure 11 
is based on that of random tumbling in free molecular flow with 
the equivalent diameter of one inch. The shape factor B has 
significant effect on total recession, while both y and geometry 
style (long or short cylinder) have only negliglble^effects . 

Figure 12 shows the same conclusion for the total recession cal- 
culations based on the ballistic coefficient on continuum flow. 


j^^caij be approxlmimately calculated by the following 


\ 
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CONCLUSION 


This study indicates that surface roughness effect on the 
total drag of entry tektite is negligible, and can be ignored 
for tektite studies, Por shapes other than sphere, the tumbling 
of tektite has significant effect on the total drag and on the 
consequent entry trajectory. In addition, the shape factor B 
was found to be a good parameter for the measurement of irregular 
shape and its tumbling effect on tektite entry. Rough estimation 
of the tumbling effect on the total recession can be estimated 
once the shape factor B and the size of tektite are given. 
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TABLE 1 


Relative Ballist Coefficient of a Cylinder with Respect 
to the Equivalent Sphere (/^sphere / /^cylinder) 


Flow Regime 


Case 

Side -on 

End -on 

Tumbling in 
pitch plane. 

Random 

tumbling 


Free Molecular 

J 1/3 
0.9717 (-) 
cl 

d 2/3 
0. 7632 (-^) 


0.6186 0 + 0.4858 Q 

d li 

nin d2/3 

0.7852 +0,3816 (^) 


Continuum 

j2 1/3 

1.2955 

d2/3 
1.3875 (;g) 

? 1/3 d 2/3 

0. 5497 (J) + 0. 5888 (j) 

Pl/3 d2/3 

0.7773 (J + 0. 3469 {|) 
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FIGURE I - VARIATION 
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FIGURE 4 - TUMBLING EFFECT AND SHAPE FACTOR IN CONTINUUM FLOWS 
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FIGURE 6 - SKIN-FRICTION DRAG COEFFICIENT AS A FUNCTION OF THE “INTERACTION 
PARAMETER (FROM REFERENCE 9) 
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TUMBLING EFFECT ON TOTAL RECESSION (CONTINUUM) 
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